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Abstract 

A mathematical model has been developed to determine heat transfer during vapor deposition of source materials under a 
variety of orientations relative to gravitational accelerations. The model demonstrates that convection can occur at total 
pressures as low as 10” 2 mm Hg. Through numerical computation, using physical material parameters of air, a series of 
time steps demonstrates the development of flow and temperature profiles during the course of vapor deposition. These 
computations show that in unit gravity vapor deposition occurs by transport through a fairly complicated circulating How 
pattern when applying heat to the bottom of the vessel with parallel orientation with respect to the gravity vector. The model 
material parameters for air predict the effect of kinematic viscosity to be of the same order as thermal diffusivity, which is 
the case for Prandtl number ~ 1 fluids. Qualitative agreement between experiment and the model indicates that 
6-(2-methyl-4-nitroanilino)-2,4-hexadiyn-l-ol (DAMNA) at these pressures indeed approximates an ideal gas at the 
experiment temperatures, and may validate the use of air physical constants. It is apparent that complicated nonuniform 
temperature distribution in the vapor could dramatically affect the homogeneity, orientation, and quality of deposited films. 
The experimental test is a qualitative comparison of film thickness using ultraviolet-visible spectroscopy on films generated 
in appropriately oriented vapor deposition cells. In the case where heating of the reaction vessel occurs from the top, 
deposition of vapor does not normally occur by convection due to a stable stratified medium. When vapor deposition occurs 
in vessels heated at the bottom, but oriented relative to the gravity vector between these two extremes, horizontal thermal 
gradients induce a complex flow pattern. In the plane parallel to the tilt axis, the flow pattern is symmetrical and opposite in 
direction from that w-here the vessel is positioned vertically. The ground-based experiments are sufficient preliminary tests of 
theory and should be of significant interest regarding vapor deposited films in microgravity. 


1. Introduction 

In recent years, a great deal of interest has been 
directed toward the use of organic materials in the 


Corresponding author. Fax: 4- 1 205 544 2102. 


development of higher efficiency optoelectronic de- 
vices. There is a myriad of possibilities among or- 
ganics which allow flexibility in the design of unique 
structures with a variety of functional group objec- 
tives. The use of nonlinear optical (NLO) organic 
materials such as thin film waveguides allows full 
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exploitation of their desirable qualities by permitting 
long interaction lengths and large power densities 
with modest power input [1]. There are several meth- 
ods in use to prepare thin films, such as Langmuir- 
Biodgett [2-4], growth from sheared solution or melt 
[5,6], and melt growth between glass plates [7], 

Epitaxial growth on ordered organic and inorganic 
substrates and variations in processing conditions 
have been useful for preparing highly oriented poly- 
diacetylene (PDA) and phthalocyanine (Pc) films 
[8-10]. The degree of significance relating process- 
ing conditions to uniformity in thickness, degree of 
orientation, and optical properties for a specific pro- 
cessing technique is the general focus of the work of 
the research group at the Marshall Space Flight 
Center (MSFC). A study on the effect of processing 
conditions relevant to thin film deposition by physi- 
cal vapor transport (PVT) is particularly difficult 
because of the possibility that convection may play a 
major role. It is a goal of some researchers having an 
interest in PVT to produce good quality anisotropic 
films, therefore, an important yet understudied re- 
quirement should be to assess the role of gravity 
during processing. This may be particularly true for 
the vapor deposition of diacetylenes where subse- 
quent polymerization in the crystal is topochemical 
and occurs readily only when neighboring monomer 
molecules are sufficiently close and suitably oriented 
[11]. Likewise, this requirement is equally viable for 
the vapor deposition of Pc’s in view of the results of 
microgravity experiments by the 3M Corporation 
involving the preparation of thin films of copper Pc 
(CuPc) [12-22]. Indeed, a variety of microstructural 
forms were obtained in thin films of CuPc dependent 
on processing methods and conditions. Small changes 
in processing parameters caused large changes in 
film orientation. Microgravity grown CuPc had sev- 
eral desirable features, which indicate that the growth 
of organic films in low-g may result in better quality 
films for optical and electrical applications [21,22]. 
The dramatic 3M result was generally unexpected 
and has been the source of considerable discussion 
over the 4-6 years since the performance of these 
experiments. 

At MSFC, we have begun studying the thin film 
growth of metal-free Pc [23], CuPc, and a polymer- 
ized 2-methyl-4-nitroaniline diacetylene derivative 
(PDAMNA) [8,24]. Experimentation on the deposi- 


tion of these materials include PVT processing, al- 
though the conditions at which transport occurs dif- 
fer considerably between the various materials. The 
deposition of the diacetylene monomer (DAMNA) 
occurred at a pressure of 10~ 2 mm Hg. Most argu- 
ments, usually of a qualitative nature, suggest that 
convection is not a factor at pressures below about 
10 mm Hg, because momentum transfer is arguably 
minimal due to large mean free paths. This would 
suggest that most processes could be designed to 
operate at low enough pressures as to negate orienta- 
tional and other possible gravitationally sensitive 
effects. Secondarily, our purpose is to begin a sys- 
tematic experimental and theoretical study to deter- 
mine a range of pressure thresholds for w^hich one 
would need to be concerned with convection. For Pc, 
the 3M method [9] involved deposition of metal-free 
Pc as substrate on the ground at about 10 -f1 mm Hg. 
At these very low pressures, orientation of the 
molecular /;-axis was only a function of temperature 
with no gravitational dependence. Epitaxial growth 
of CuPc at much higher pressures onto the metal-free 
substrate required the low gravity environment for 
maximum orientation of the CuPc film (PVT for 
organic solids (PVTOS)). 

A major goal of our research program is to under- 
stand factors for improving film quality and optical 
properties in vapor deposited films. Important as- 
pects of any study involving fluids, as in vapor 
transport, are driving mechanisms for heat transfer 
with natural convection and diffusion processes 
which determine flow profiles and temperature dis- 
tributions. In this study, we have considered initia- 
tion and development of time-dependent, three-di- 
mensional fluctuations of flow profiles and tempera- 
ture distribution. These phenomena are driven by 
heat transfer in parallel and oblique directions wdth 
respect to gravitational acceleration. Furthermore, 
consideration of pressure effects suggest that in unit 
gravity, convection can occur at lower pressures than 
previously thought. The ground-based portion of this 
study investigates factors influencing vapor transport 
varying reactor orientation with respect to the gravity 
vector. We seek to correlate cell orientation and 
other processing parameters with film quality. Exten- 
sive experimental and theoretical work will provide 
data to assess the value of future microgravity pro- 
cessing. The primary objective of this work, how- 
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ever, is to concentrate on the theoretical framework 
which describes heat flow and thermal profiles using 
time-dependent, three-dimensional formulations at 
various cell orientations. Tests of theoretical predic- 
tions initially rely on relative film thickness mea- 
surements following depositions in appropriately po- 
sitioned cells. 

Recognizing the criticality of selected boundary 
conditions, the size and geometry of the container in 
these simulations represent initial parameters and 
define actual but roughly designed ground-based ex- 
periments. The results of this study will allow refine- 
ment of cell dimensions and design for assessment of 
past and future ground-based work, and provide an 
additional basis for investigating microgravity pro- 
cessing by PVT. 

2. Mathematical formulation 

In this section, we introduce the governing equa- 
tions, initial and boundary conditions to assess the 
flow fields and temperature distributions upon appli- 
cation of a thermal gradient to an ideal gas. There is 
no species diffusion component to the formulation 
that accounts for boundary fluxes driven by sublima- 


z 



Fig. 1 . Geometrical illustration of coaxial cylinders to he used for 


tion at the source and condensation at the sink 
[24-26]. Subsequent computations will contain 
species diffusion equations with inclusion of advec- 
tive-diffusion flow. These computations approach 
steady state and for unit gravity, time-dependent 
flows and temperature distributions generally result 
from initial transients. 

We did not seek to observe oscillatory flows. 
Time-dependent computations during microgravity 
processing, however, can serve the important pur- 
pose of correlating experimental results with fluctu- 
ating accelerometer data. Subsequent computations 
will explore time-dependent oscillatory flow behav- 
ior with Fourier component analyses for both unit 
and microgravity environments. 

Consider a coaxial cylinder with an inner radius 
of /?, and an outer radius of /? 2 . The heights of the 
outer and inner cylinders are L 2 , and (L 2 — L,), 
respectively. Fig. 1 shows the geometrical configura- 
tion of the coaxial cylinder. For the purpose of 
considering the mathematical formulation, we adopt 
cylindrical coordinates (t\6,z) with corresponding 
velocity components and corresponding 

components of gravity acceleration (g r ,g H ,g : ). Figs, 
la and lb show the geometrical configuration of a 
coaxial cylinder in an r-z plane and the r-0 plane. 



(b) 

mathematical modeling, (a) an r-z plane, and (b) the r-0 plane. 
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respectively. There are two temperature zones in this 
coaxial cylinder. During our growth of DAMNA 
films [27], the external walls of the outer cylinder are 
surrounded by an oil jacket, which maintains a con- 
stant temperature, ~ 120°C, while the walls of the 
inner cylinder keep a constant temperature of ~ 
30°C. The cell, containing air, is evacuated to a 
pressure of 10 2 mm Hg (= 1.333 N/m 2 ). The air 
fills the gap between the outer and inner walls of the 
coaxial cylinder. The How profiles and temperature 
distribution of the air in the cell at relatively low 
pressure are influenced by the natural convection of 
the flow pattern in this problem. The governing 
equations can be shown as follows: 

Continuity equation, 


1 a 

r dr 


1 dr chv 

( ru) + + — = 0. 

v 7 r d0 dz 


(i) 


Momentum equations, 

/ du du 

+ U + 

\ dt dr 

dp 
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dw \ 
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(3) 


(4) 


and 


U, 

— g(sin ip cos 0 , — sin ip sin 0, —cos ip ) . 


Energy equation. 


P c p 


a t a t r a t ar\ 

dt ar r dO dz J 


kV 2 T. 


(5) 


Here, p, p, 7\ p,, k , and g denote density, pressure, 
temperature, dynamic viscosity, thermal conductivity 
of air, and background gravity acceleration, respec- 
tively. ip is the angle between the gravity vector and 
the axis of the cylinder. In this case, and p x 
represent the air temperature and density at 10“ ~ 
mm Hg and heat sink at the wall of the inner 
cylinder at T = 30°C. In these mathematical expres- 
sions, the Boussinesq approximation is employed to 
formulate the relationship between local fluctuations 
of density and temperature through the thermal ex- 
pansion coefficient which replaces the equation of 
state. 


3. Initial and boundary conditions 

The initial and boundary conditions have been 
assigned explicitly to solve governing Eqs. (t) to (5) 
simultaneously. The initial condition is at t = 0, all 
the flow profiles vanish, i.e., 

u = r = vv = 0 everywhere. (6) 

The temperature of fluid in the inner cylinder is 
30°C. The temperature of fluid in the water jacket 
surrounding the outer cylinder is 120°C. These con- 
ditions are expressed by, 

T = 30°C at r</?,,0< 0<2 tt, L x <z<L 2 , (7) 


T= 1 20°C at z = 0\ r>0;0< <2 tt, (8) 

at z = L 2 ; r > 0 ; 0 < 0<2 t7, (9) 

at 0 < z < L 2 ; r = R 2 ; 0 < 0 < 2tt . (10) 

The pressure profile inside the container at t — 0 is 
p= 10“ 2 mm Hg = 1.333 N/m 2 . (11) 


The boundary conditions are that at time / > 0, 
the temperature distribution shown in Eqs. (7)-(10), 
and pressure profiles shown in Eq. (II) remain the 
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same. For flow profiles, no penetration and no slip 
conditions along the surfaces of the solid wall shall 


apply, namely, 

u = i=w = 0 (12) 

along the surfaces of 

r = /? I ,O<0<27T, L, <z<L 2 , (13) 

r > 0, 0 < 6 < 2tt, z = 0, (14) 

r > 0, 0 < 0 < 2 t7\ z = L 2 , (15) 

r = /? 2 , 0 < 6 <2tt i 0 < z < L 2 . (16) 


The coaxial cylindrical container, in Fig. 1, has the 
following specifications: L { = 6 cm, L 2 = 12 cm, 
R { = 0.9 cm, and R 2 = 1.5 cm. The physical param- 
eters for determining flow fields and temperature 
distribution were gravity acceleration g 0 = 9.81 
m/s 2 , air pressure = 1 .333 N/m 2 , air density = 
1.71 X 10 5 kg/m 3 , air kinematic viscosity coeffi- 
cient = 1.17 m 2 /s, coefficient of thermal- expansion 
for air p = 2.87 X 10~ 3 1/K and Prandtl number for 
air = 0.72, heat conduction coefficient for air = 1.8 
X 10 " 2 J/(m s K), and the constant pressure spe- 
cific heat, o p , for air = 1.014 X 10 3 J/(kg K). 

Based on a Boussinesq approximation, application 
of the thermal expansion formulation is appropriate 
due to the considerable temperature and density 
changes in the ampoule that arise from significant 
non-isothermal conditions. The most significant driv- 
ing parameters in this study, avoiding an exhaustive 
discussion of dimensionless parameters, are the iner- 
tia, buoyancy and viscous terms, 

ww/L - g z (SAT ~ v w/df , (17) 

respectively, which are on the same order of magni- 
tude in the momentum equation shown in Eq. (4); 
and convection and conduction terms, 

w/LAT ~ aAT/8 ~ , (18) 

respectively, which are also on the same order of 
magnitude in the energy equation shown in Eq. (5). 
Here, 8 T is the thermal boundary layer thickness; L 
is the height of the container along the z-axis; v~ 
fx/p is the kinematic viscosity; and a = k/pc p is 
the thermal diffusivity. From convection and conduc- 
tion terms in Eq. (18); from inertia and buoyancy 
terms in Eq. (17); and from buoyancy and viscous 


terms in Eq. (17), we can obtain convective velocity 


of flows as follows: 

w-aL/fif, (19) 

w~{g z piAT)' /2 , (20) 

w~(g : 08ilLT)/ V . (21) 

By equating Eqs. (19) and (21), we have the relation- 
ship 

<5 T /L ~ Ra - 1/4 (22) 

a result of the balance between buoyancy and vis- 

cous forces. By equating Eqs. (19) and (20), we have 
another relationship 

8 t /L~ ( Ra Pr) 1 4 - (Gr) l/4 (Pr) 1/2 (23) 


based on the balance between buoyancy and inertia 
forces. Here the Rayleigh number (Ra), Grashof 
number (Gr), and Prandtl number (Pr) are defined as 


Ra = g z pL^AT/av, (24) 

Gr = g. pL^AT/v 2 - Ra/Pr, (25) 

and 

Pr ~ v/ tx . (26) 

Since the heat transfer coefficient scales as k/8 T 
[25], the Nusselt number (Nu) varies as 

Nu = hl/k - Ra 1/4 (27) 

from Eq. (22) based on the balance between buoy- 
ancy and viscous forces inside the thermal layer; and 
Nu varies as 

Nu - (Ra Pr) 1/4 ~ (Gr)' /4 (Pr)' /: (28) 


from Eq. (23) based on the balance between inertia 
and buoyancy forces inside the thermal layer. As 
indicated by Bejan [28] both theoretical and labora- 
tory measurements show that Eq. (27) is particularly 
applicable for Pr > 1 fluids, while Eq. (28) is domi- 
nated by Pr I fluids. For the case of the present 
research, the dominant parameter is the Rayleigh 
number rather than Grashof number [28] since for 
air, Pr ^ 1 . 

In this study, it is shown in the numerical simula- 
tion that for the range of parameters considered, the 
maximum induced flow velocity driven by convec- 
tion is 0.165 m/s. The Reynolds number (Re = 
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wD/v, where D is the diameter of the container) for dent, three-dimensional flow profiles and tempera- 

the flow pattern based on this flow velocity is 0.025 ture distributions numerically. The computer algo- 

m/s, which is far below the threshold for production rithms employed in this study have been developed 

of turbulent flow fields. The Grashof and Rayleigh and illustrated in our earlier studies [29-39] and will 

numbers for the present case study are 0.0184 and not be repeated in this paper. The time step is 

0.0133, respectively, for unit gravity. The aspect determined automatically based on the size of grid 

ratio of the reaction chamber (height : cross-sectional points and the velocity of flow fields. Since the 

radius = 6: 1.5) is sufficiently large such that any thickness of the boundary layer is inversely propor- 

radial gradient (T - Tj will result in buoyancy- tional to the square root of Reynolds number, the 

driven flows [28]. That is, there is no critical Rayleigh spacing is adjusted such that it decreases with in- 
number requirement for the onset of convection in creasing Reynolds number [40,41]. 

the present study. In this study, two sets of grids, with grid sizes 

17 X 34 X 42 and 23 X 34 X 82 along (r,0,z) coor- 
dinates, have been designed to carry out numerical 
4. Mathematical simulation of flow profiles and computations. For a Reynolds number of 0.025 with 
temperature distributions maximum induced velocity of 0.165 m/s, as in the 

present case, there is concluded to be no difference 
Eqs. (l)-(5) subject to initial conditions shown in for numerical results obtained for either grid chosen. 

Eqs. (8)— ( 1 1 ), and boundary conditions shown in By using the grids 17 X 34X42 in the computation 

Eqs. (12)-(16), have been used to solve time-depen- of flows with higher Reynolds number (such as 


(a) t » 3.3 s (b) t = 20 s (c) t = 30 s (d) t - A0 s 

m 4 Sfi nm/q = ft fil rm/c \/~ = ft 7ft /« V™ . ft 7R rm/s 



Radius (cm) Radius (cm) Radius (cm) Radius (cm) 

Fig. 2. Flow profile initiation and development along the c-axis parallel to the plane 0 = 0°-180° and to gravitational acceleration (ip = 0°). 
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Re > 1500) and induced velocity of 1.66 m/s and 
higher, extremely severe numerical instabilities of 
the flow field result. With grids 17 X 34X42, the 
average time step determined from the size of grid 
points and velocity of flow fields is 0.00076 s, while 
the CPU time required by the Supercomputer CRAY 
II to execute one time step is 3.019 s. The total CPU 
time consumed for each of the three cases in this 
study is 12 000 s. The cell orientation varies relative 
to gravity according to 1^=0°, 45°, and 180°, by 
rotating the z-axis. We observe flow profiles in the 
cell in the orthogonal planes, 0 = O°-18O° and 0 — 
90°-270° to investigate the characteristics of convec- 
tion-driven flows and temperature profiles. Various 
references have been cited by Bejan [28] in the study 
of induced convective heat transfer with fluid layers 
heated from below, top, side, oblique, etc., under 
normal atmospheric pressure. 

As mentioned in the Introduction, this study in- 
vestigates the time dependence from initiation and 
development of flow fields and temperature distribu- 
tion throughout the vapor deposition process. The 
objective is to determine the sensitivity of the forma- 
tion of NLO thin films to convection in the relatively 
low pressure environment of 10“ 2 mm Hg. To our 


knowledge, there are no studies of induced convec- 
tive heat transfer at relatively low pressure, hence, 
there is no available literature for comparison. How- 
ever, the following discussion concludes there is 
qualitative agreement on flow profiles and tempera- 
ture distribution for induced convection at both nor- 
mal atmospheric pressure and 10 2 mm Hg. 

5. Convection driven by processing with the cell 
oriented parallel ( if/ = 0°) to gravitational acceler- 
ation 

Fig. 2 shows the time dependency for the devel- 
opment of flow fields (V m = maximum velocity at 
time, t) upon bottom heating of the cell oriented in 
the direction parallel (tp — 0°) to gravitational accel- 
eration. Fig. 3 shows time animation of the initiation 
and formation of temperature profiles due to induced 
convection in the r-z plane, parallel to the plane 
61 = 0°- 180°. The curves shown in Fig. 5 are the 
isothermal lines with temperatures in °C at the corre- 
sponding locations. Figs. 2 and 3 can be summarized 
as follows: (a) At time t = 0 (not shown), there is 
very sharp temperature variance from 30 to 120°C in 



Fig. 3. Initiation and formation of the temperature profile along the c-axis parallel to the plane 6 = 0°- 180° and to gravitational acceleration 
( ip = 0°). 
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(a) t - 3.3 s 
V m = 4.56 cm/s 



(b) t = 40 s 
V m = 8.78 cm/s 



Radius (cm) 


Radius (cm) 


Fig. 4. Flow profile initiation and development along the “-axis 
parallel to the plane 0 = 0°-l80 c with top healing (i// = 180°). 


a fairly thin layer, (b) An initiation of downward 
flow from the bottom of the central column toward 
the bottom of the outer cylinder immediately creates 
a gradually expanding positive temperature gradient 
along the downward directions for time i > 1.5 s. (c) 
A pattern of induced inflow from the surface of the 
outer cylinder toward the center promotes a parabolic 
temperature profile with the downside tip point 
around the center. This pattern is due to downward 
cold streams from the bottom of the central column 
for time t > 1.5 s. (d) An induced upward profile of 
temperature distribution around the surface of the 
outer cylinder is a result of an upward flow pattern 
along the corresponding surface, (e) A downward 
elongation profile of temperature distribution imme- 
diately around the neighborhood of the surface of the 
inner cylinder (central column). This pattern is a 
result of the induced downward flow along the corre- 
sponding surface. 


6. Convection driven by processing with the cell 
oriented in the reverse parallel direction (i jf — 
180°) to gravitational acceleration 

Fig. 4 and Fig. 5 show the time dependency for 
the development of induced flow fields upon top 
heating of the cell oriented in the parallel direction 
(t/>= 180°) to gravitational acceleration. This prob- 
lem can be divided into two sections. For z > 6 cm, 
the problem can be associated with horizontal heat- 
ing from the outer walls of the cylinder toward the 
central wall within a narrow cylindrical slit of space. 
For z < 6 cm, the problem is that of vertical heating 
from the top (downward gravity direction). As the 
source material is located at the surface of z = 0 
(now- at the top), and the location of vapor deposition 
is at the bottom horizontal surface of the central 
column at z. = 6 cm, the flow profile at c > 6 cm has 
very little effect on the flow pattern at the deposition 
surface. Thus, we can ignore the flow profiles in- 
duced within a narrow cylindrical slit of space. As to 
the flow at z < 6 cm. Fig. 4 shows a stable stratified 




Fig. 5. Initiation and formation of the temperature profile along 
the :-axis parallel to the plane 0 = O°-18O° with top heating 
(<// = 180 °). 
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medium and natural convection is absent provided 
there is no side heating. 

Fig. 5 shows the time evolution of initiation and 
formation of temperature profiles for this case. The 
curves in Fig. 5 are the isothermal lines with temper- 
ature in °C at the corresponding locations. Again, the 
figure clearly shows two regions. For the region 
z > 6 cm, horizontal heating from the side walls 
drive the convective flow toward the central wall 
within a narrow cylindrical slit. Fig. 5 depicts a 
series of isothermal curves with decreasing tempera- 
ture from the side wall toward the central column. 
As indicated earlier, this region is not important in 
our study of vapor deposition from source materials 
located on the surface at z = 0 and transferring to 
the bottom of the central column at z = 6 cm. There 
is not much variation in the dynamics of temperature 
distribution for the region z < 6 cm because the 
region is filled with the stable stratified medium, and 
natural convection is completely prohibited. Never- 
theless, the heating and consequential flow transient 
is of interest. 


(a) t - 3.3 s 



(b) t - 40 3 
V m = 8.78 cm/s 
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7. Convection driven by the oblique direction at 
45° to gravitational acceleration 

In this case, the cell z-axis is allowed to tilt at 45° 
relative to gravitational acceleration parallel to the 
plane, 0 = O°-18O°. Fig. 6 illustrates convective 
flows induced by a vector summation of gravita- 
tional acceleration with a horizontal component at 
0.707g 0 (g D sin 45°) rightward, and vertical compo- 
nent at O.707g o (g 0 cos 45°) downward. 

Comparison of Fig. 2 and Fig. 6, and Fig. 2 and 
Fig. 7 show the differences in induced flow fields in 
the 9 = 0°- 1 80° and 9 = 90°-270° planes, respec- 
tively, driven by vertical and titled cavities. There 
are two recirculation flows generated with centers at 
(r c ,z c ) = (± 1.90, 2.7) cm and at (r t .,z c ) = ( ±0.8, 
3.2) cm for flows induced by vertical and tilted 
orientations, respectively. For the vertical orienta- 
tion, the flow depicted in Fig. 2 is axisymmetric, 
hence also representative of flows in the 9 = 90° - 
270° plane. The flows are upward along the surface 
of the outer cylinder and downward from the central 
column. These are Benard-type flows driven by 
counter-rotating cells (clockwise on the left and 


Radius (cm) Radius (cm) 

Fig. 6. Flow profile initiation and development along the ->axis 
parallel to the plane 0 = 0°- 180° and oblique to gravitational 
acceleration (t/r = 45°). 

counter-clockwise on the right) due to an incipient 
instability in the narrow cylindrical cell cavity. There 
is no critical Rayleigh number for convection, which 
assumes an infinite extent of the cell width. For the 
tilted cell, a resultant asymmetric flow in the 9 = 0°- 
180° plane approaches an antisymmetric flow profile 
resulting from differential heating between vertical 
and opposite walls of a cavity. That is, an incipient 
instability causes heat to flow along the bottom 
surface (now tilted upward), and upward along the 
side wall, and along the cold surface (downward), 
and downward along the other side wall. The two 
recirculation flows for this orientation appear in the 
0=9O°-27O° plane (Fig. 7). There is an upward 
flow from the bottom surface of the outer cylinder 
transporting heat toward the bottom surface of the 
central column with downward flows on the outside 
walls (clockwise on the right and counter-clockwise 
on the left). This asymmetric three-dimensional flow 
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profile is quite complex and represents a greater 
degree of convection in the cell cavity for the tilted 
orientation. 

Fig. 8 and Fig. 9 show time evolution of the 
initiation and formation of temperature profiles along 
the z-axis parallel to 0 = 0°-180° and 0 = 90°-270°, 
respectively. Fig. 8 depicts the temporal variation of 
isothermal lines as a result of the tilted orientation 
with respect to gravity. The isotherms are consistent 
with the motion of flow fields. Comparison of Fig. 6 
for the velocity profile and Fig. 8 for the temperature 
profile demonstrates that the flow field affects the 
temperature field. Fig. 9 shows similar temporal 
variation of isothermal lines from another point of 
view. Initiation and development of a cone-shaped 
temperature profile is the result of direct upward 
flow from the bottom surface of the outer cylinder 
transporting heat toward the bottom surface of the 
central column. This is viewed in the 0 = 9O°-27O° 


plane and is consistent with the flow profile in Fig. 

7. These preliminary observations suggest that vari- 
ous convection patterns can affect film thickness and 
quality during deposition by vapor transport. 

8. Near ultraviolet / visible absorbance spectra of 
vapor deposited diacetylene-MNA (DAMNA) films 
as a function of cell orientation 

In this work, it is helpful to utilize Beer-Lambert’s 
relationship for transmission of radiation through a 
medium to test gravitationally sensitive flow pattern 
predictions. We would expect that the flow pattern in 
Fig. 7 (cell tilted at 45° relative to gravity vector) 
would affect film quality over most of the deposition 
surface differently than the flow pattern depicted in 
Fig. 2 (cell oriented parallel to gravity vector). Pre- 
liminary experiment suggests this to be the case with 


(a) t « 3.3 s 
V m = 4.56 cm/s 



(b) t = 20 s 
V m = 8.61 cm/s 
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(c) t * 30 s 
V m = 8.78 cm/s 



(d) t =* 40 s 
V m « 8.78 cm/s 



Radius (cm) 


Radius (cm) 


Radius (cm) 


Radius (cm) 


Fig. 7. Flow profile initiation and development along the c-axis parallel to the plane 6 — 90°— 270° and oblique to gravitational acceleration 
( ip ~ 45°). 







298 


D.O. Frazier et ai /Journal of Crystal Growth 171 (1997) 288-302 


respect to film thickness. The Beer-Lambert law 
relative to film thickness can be written 

/ = / 0 e-“ / , (29) 

where / 0 is incident radiation intensity, I is the 
reduced radiation intensity after passing through the 
film, a a proportionality constant containing the 
molar absorption coefficient, e, and concentration of 
absorber, and / is film thickness. The molar absorp- 
tion coefficient represents a molar cross section for 
absorption. The greater the cross section of the 
molecule for absorption and absorber concentration, 
the greater the attenuation of the intensity of the 
beam. Likewise, film thickness, / also attenuates 
the beam intensity accordingly. The ratio /// 0 is a 
measure of beam transmittance, 7. We would expect 
a larger intensity in the absorption bands of films 
deposited in cells tilted at 45° over those where 
deposition occurred in cells positioned vertically. We 
may define the dimensionless product A — — a / as 
absorbance, which incorporates all of the contribu- 
tors to beam attenuation. Since the nature of the 


(a) t = 3.3 s (b) t = 20 s 



Fig. 8. Initiation and formation of temperature profile along the 
acceleration (tj/= 45°). 


material, DAMN A, and concentration (pure material) 
are identical in both cases, then absorbance is only a 
function of film thickness. Fig. 10a is a wavelength 
scan in the range 190 to 820 nm of a vertically 
deposited film, and Fig. 10b is that of an obliquely 
deposited film. Table 1 contains representative ab- 
sorbance intensities from scanning an approximately 
1 cm diameter spot in similar vicinities of films 
formed during vapor deposition in vertically and 
obliquely oriented cells. Visually, the films from the 
obliquely oriented cells are a deeper yellow color 
than those from vertically oriented cells. We qualita- 
tively assume generally thicker films from the rela- 
tive appearances. The surfaces of these monomer 
films are translucent and microcrystalline. At the 
present time, it is also preferable for us to discuss 
relative film thickness from beam attenuation in 
qualitative terms (although more quantitative than 
visual observation) while observing that the spectro- 
scopic irradiation spot sizes are large enough to 
average over about 50% of the film surfaces. The 
same relative result occurred, in that beam attenua- 


(c) t = 30 s (d) t = 40 s 



c-axis parallel to the plane # = 0°-180 o and oblique to gravitational 
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tion was greater at each maximum absorbance wave- 
length repeatedly. We may consider Table 1 a quali- 
tative representation of a consistent result. Deposi- 
tions were of 30 min duration. Cell dimensions were 
the same as those discussed in the modeling section. 
The cells were pumped to 10 2 mm pressure prior to 
deposition, and sealed. Without the presence of 
DAMNA in the vessel, appreciably lower total pres- 
sures are obtainable. We may infer that the evacua- 
tion pressure reaches a minimum because of the 
vapor pressure of DAMNA, which will be equal to, 
or greater than, the evacuation pressure, 10 -2 mm. 
From a different study, the measured vapor pressure 
of 4-N,N-dimethylamino-4'-nitrostilbene (DANS) at 
120°C was reportedly 0.374 mm [42]. This is a 
relatively large organic molecule having a molecular 
weight of 21 1 g/mol as compared to DAMNA with 
a molecular weight of 247 g/mol. Considering struc- 
tural and size similarities, we may approximate simi- 
larities in vapor pressure. Indeed, heptadecanol 
(molecular weight = 256.5 g/mol), also an alcohol 
with possible hydrogen bonding in its condensed 
phases such as expected of DAMNA, has a vapor 
pressure of 10” 1 mm at 120°C [43]. There is no 
reason to expect DAMNA to differ drastically from 
these measured vapor pressures, and we may approx- 
imate that the vessel evacuation pressure minimum is 
largely due to the vapor pressure of DAMNA. 

The vacuum was virtually unchanged during de- 
position. The oil baths circulated at 120 and 30°C 
(source and sink, respectively) prior to allowing flow 
into the experimental cell. The temperature differen- 
tials between the baths and respective oil jackets of 
the experimental cells were small ( ~ 2°C). This 
allowed operation during deposition at constant heat 
fluxes of ~ 120°C at the source and ^ 30°C at the 
sink, consistent with theoretical conditions. Since the 
mathematical formulations only considered the ther- 
mal flow profiles using air data, these results are 
estimates. That is, the effects of mass transport are 
not included in these assessments. However, at low 
pressures and moderately high temperatures, we can 
expect ideal gas behavior as a first approximation; 
namely, that gaseous molecules are essentially point 
masses. In comparing experiment with the model, we 
expect that heats of sublimation and condensation 
will have little effect on the thermal profile. Vapor 
phase estimates on the basis of ideal gas considera- 


<a) t = 3.3 s (b) t = 40 s 



Fig. 9. Initiation and formation of temperature profile along the 
--axis parallel to the plane 0 = 90°-270 c and oblique to gravita- 
tional acceleration (t // = 45°). 

tions and DAMNA condensed volume, places 
DAMNA’ s mole fraction in the available cell cavity 
to be very high relative to air (^ 1). The model 
material parameters for air predict the effect of kine- 
matic viscosity to be of the same order as thermal 
diffusi vity , which is the case for Pr ~ 1 fluids (Eq. 
(26)). The classical heat capacity for a gaseous 
molecule, such as DAMNA, with its large number of 
vibrational degrees of freedom, is only achieved at 
high temperatures. That is, the heat capacity of 
DAMNA will approach that of a diatomic molecule, 
such as the major components of air, since its vibra- 
tional modes are unable to store energy at these low 
temperatures. The thermal diffusi vity approximation 
using air data is, therefore, a fair one assuming 


Tabic 1 

Representative absorbance intensities of deposited films 


A (nm) 

A m;ix (cell vertical film) 

(cell oblique film) 

422 

0.200 

0.385 

214 

0.395 

0.546 
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similar heat conduction coefficients. In the classical 
limit, the higher heat capacity would yield Pr » 1 . 
Because an ideal heat capacity approximation is an 
overestimate at the operating temperatures, deviation 
from ideal behavior causes Pr ^ 1 , closer to that of 
air. Without actual data for DAMNA, we are left 
with experimental data from the low operating pres- 
sures and temperatures to compare with modeling 


estimates. Close agreement between experiment and 
the model would indicate that DAMNA at these 
pressures indeed approximates an ideal gas with the 
appropriate deviations which tend toward validation 
of the use of air physical constants. Furthermore, the 
constant heat flux provided by the circulating bath is 
a factor in neutralizing thermal effects specific to 
individual molecules. In keeping with the assump- 




Fig. 10. Representative ultraviolet/visible spectrum of a deposited DAMNA film in (a) the deposition cell oriented vertically with respect to 
the gravity vector, and in (b) the deposition cell oriented at 45° with respect to the gravity vector. 
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tions, the mathematical model correctly predicts 
qualitative differences in film properties between 
vertically and obliquely oriented cells. That is, more 
convection in cells having the tilted orientation is 
apparently responsible for correspondingly greater 
film thickness. Future work will strive to correlate 
experimental kinetics of film growth with calculated 
heat and mass transfer rates. 


9. Discussion and conclusion 

The purpose of the present study is to investigate 
heat transfer characteristics in an ideal gas in a 
thermal gradient at relatively low pressure and vari- 
ous orientations with respect to the gravity vector. 
Mathematical formulations yield fluid profiles, air 
and temperature distributions driven by heat transfer 
processes during natural convection at orientations 
parallel and oblique with respect to gravity. Several 
experiments to grow films by PVT have resulted in 
some controversy. The argument raises the question 
of minimum pressure levels for convection to be a 
factor. Neither has it been clear as to the types of 
convection which might occur during PVT that would 
affect the quality, homogeneity and orientation of 
vapor deposited films as in the 3M microgravity 
experiment. Mechanisms thus described to discuss 
the 3M experiment are generally in terms of thermal 
stress and creep due to high horizontal thermal gradi- 
ents at the corners and thermal gradients along the 
walls [44]. Although our ground-based experiment 
only roughly approximates PVTOS, Fig. 8 illustrates 
the presence of such gradients in an r-z plane 
parallel to the plane, 0 = 0°-180°. 

In this study, numerical computation of a full set 
of nonlinear fluid dynamical equations show that the 
vapor deposition process has to be accomplished 
through a fairly complicated set of flow patterns 
(particularly for three-dimensional flows at i/y = 45°) 
of recirculation in unit gravity al pressures as low as 
10“ 2 mm Hg. Since these equations also demon- 
strate an enhancement of convection as we move 
from vertical to oblique orientation of the cell with 
respect to the gravity vector, it is feasible to apply a 
simple preliminary test of these calculations by 
macroscopic assessments of film quality of vapor 
deposited films in respectively oriented cells. The 


preliminary experimental result tends to support the 
numerical computations. 
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